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Arc Measurement of Some Argon Transition Probabilities 
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Transition probabilities and line shape parameters for three argon lines have been measured in a 
wall-stablilized argon arc containing a trace of hydrogen. The H# line shape measurement and the 
assumption of local thermodynamic equilibrium provided the determination of the argon level popu- 
lations. From these three lines ten other argon line transition probabilities have been obtained by 
relative measurements in a pure argon arc. The results are compared with those of other investiga- 
tors and sources of error are discussed. 
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1. Introduction 

The importance of argon transition probabilities 
for plasma diagnostics is reflected in the large amount 
of recent experimental work on the measurement of 
these quantities [1, 2, 3, 4]. 1 The most fruitful source 
of these measurements has been the high current 
thermal arc, but unfortunately even here in similar 
sources the results of different investigators may differ 
by factors of two, well outside the error estimates 
attributed to the measurements by their authors. 
Consequently an additional attempt to measure argon 
transition probabilities in a thermal arc source was 
undertaken. The particular technique employed 
was one which has previously been used in these 
laboratories [5]. It consists of operating a wall-sta- 
bilized arc with a gas mixture containing a trace of 
hydrogen in the gas whose transition probability is to 
be measured, and then deducing the arc plasma tem- 
perature and composition from the spectral line shape 
of the 4861 A hydrogen Balmer (3 line. From the 
plasma composition and absolute line intensity meas- 
urements the transition probabilities can be derived. 

Absolute transition probabilities were determined 
in this manner for the two AI lines at 4300 and 6965 A 
and for the All line at 4806 A. These lines were 
selected for their prominence, isolation from their 
neighbors, symmetry, relative freedom from self ab- 
sorption, and location in wavelength regions favorable 
for photoelectric detection. Some additional tran- 
sition probabilities were measured relative to one or 
another of these three reference lines in pure argon 
arcs. However, measurements on these lines are less 
satisfactory either because the lines are weak or be- 
cause of ambiguity in the location of the background. 
Most of the prominent blue argon lines in this source 
are afflicted with this background problem due to the 
closeness of strong neighbors or asymmetric line 



shapes (see fig. 1) either of which can result in an 
uncertainty in the total intensity to be ascribed to the 
lines. 

All of the data recording, handling, and calculations 
including the graphical presentation of the results 
were done by automatic machine methods to eliminate 
as many sources of error as possible. With the statis- 
tically large samples made possible by these methods 
along with elimination of subjective measurements 
and increased precision, it becomes possible to at- 
tempt a more thorough analysis of the errors and un- 
certainties involved in these measurements. 
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FIGURE 1. Atmospheric pressure argon arc spectra after Abel in- 
version [ref. 19]. 

The total intensity of AI 4300 A yields an electron density of 6.9 X 10 Ui em 1 and tem- 
perature of 12 060 °K assuming the transition probability value of table 3. 
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2. Experimental Details 
2.1. Light Source 

The radiation source is a wall-stabilized or cascade 
d-c arc of about 10 cm length and 5 mm diam. The 
arc chamber is similar to that described in reference 
[6] with the exception of improved water cooling for 
higher current operation. The graphite electrodes 
are surrounded by a separate flow of argon, preventing 
rapid electrode erosion and contamination of the test 
gas in the observation section. After passing through 
the arc the gases are exhausted to the atmosphere 
through large diameter ports. The arc character- 
istically operates stably for long periods of time for 
most test gases. For the experiments reported here 
arc currents of 40, 60, and 90 A were used. 

The arc chamber is mounted on a carriage movable 
by means of a lead screw in a direction transverse 
to the optic axis of the monochromator. An elec- 
trically actuated single-revolution clutch couples an 
electric motor to the leadscrew drive providing for 
either continuous or incremental movement of the 
arc carriage. Variation in rate of scan or increment 
size is accomplished by a set of quick-change gears 
in the drive train. 

Electrical power is supplied to the arc by either 
a bank of 100 6 V storage batteries, or a rotary con- 
verter of 500 kVA capacity and 650 V output as shown 
schematically in figure 2. Both power sources were 
used during this series of experiments with no observ- 
able difference in results. The voltage drop across 
the arc is typically 102 V in argon at 60 A, and 112 V 
with 1 percent hydrogen added to the test gas. The 
bulk of the power is dissipated in a large variable 
water resistor [7] in series with the arc supply. 

The current is automatically regulated by a solid- 
state series current regulator (see fig. 2) consisting 
of 30 power transistors in parallel, with the current 
error signal being derived from a 0.2 ft oil-cooled 
manganin shunt. The arc current is initially set and 
later monitored with a potentiometer and galvanometer 
by which variations in total arc current as small as 
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FIGURE 2. Arc circuit schematic including some details of the 
transistorized automatic current regulator. 



1 mA may be observed. A carbon pile rheostat is 
used to provide vernier adjustment of the ballast 
resistance. 

2.2. Optical and Spectrographs Arrangement 

The arc was viewed in a plane perpendicular to 
the arc axis with the image focused slightly enlarged 
on the monochromator entrance slit. A rectangular 
diaphragm opening at the lens formed the aperture 
stop of the system and limited the solid angle from the 
arc to 0.002 sr. The dispersing element was an 0.8 m 
Ebert-Fastie plane grating monochromator having a 
nominal dispersion of 10 A/mm in the first order. 
Entrance and exit slits of 15 jjl width were used in 
most of the experiments. The entrance slit height 
was masked to 1.2 mm to eliminate reflections from 
the copper plates adjacent to the observation section 
of the arc chamber and to minimize the effect of any 
slight variations of the arc column along its length. 
During the course of the experiments, RCA 931 A 
(S-4 response), Hamamatsu R-136 and RCA 7265 
(S-20 response), and ITT FW-118 (S-l response) 
photomultiplier tubes were employed to detect the 
spectral intensity. 

The monochromator grating drive is fitted with a 
single channel paper tape reader attached directly 
to the wavelength counter to allow direct digital con- 
trol of the wavelength setting. An electrical signal 
causes the monochromator to scan in wavelength 
until a punched hole in the paper tape is sensed by the 
reader. A program paper tape may therefore be 
punched to direct the monochromator through a pre- 
scribed, not necessarily evenly spaced sequence of 
wavelengths in the spectrum. The resolution of the 
tape reader is 0.05 A and the precision of the wave- 
length setting better than 0.01 A. 



2.3. Data Recording and Control System 

The large quantity of data encountered in spectro- 
metric studies involving spacial resolution of plasma 
characteristics invites the use of automatic digital 
techniques for data collection and reduction [8]. To 
facilitate data handling the response of the mono- 
chromator was recorded digitally on punched paper 
tape in a form compatible with electronic data proces- 
sing equipment. In addition, the monochromator 
wavelength drive and arc carriage positioning mech- 
anism were automatically actuated and coordinated 
so as to perform the proper variable changes during an 
experiment in the minimum time. 

In essence, the data recorder consists of amplifying, 
digitizing and tape punching equipment. The small 
current output of the photomultiplier is amplified and 
converted to a voltage signal of the order of 1 V which 
in turn is converted to a pulse train whose frequency 
is proportional to this voltage. An electronic counter 
sums the pulses for a precise interval of time, resulting 
in a digitized integral of the applied voltage over the 
counting interval. Upon command, six digits of in- 
formation from the counter are loaded simultaneously 
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into a shift register from which they may be shifted 
out serially into a code-conversion matrix and punched 
on paper tape in binary-coded-decimal form. 

Movement of the arc carriage, changes in the mono- 
chromator wavelength and recording of photomultiplier 
current are all directed and coordinated automatically 
during a run by a central programmer. The heart of 
the programmer is a 25-step, 11 level stepping switch, 
wired in a closed-loop mode; i.e., a signal from the 
preceding operation verifying its successful completion 
is necessary to initiate a succeeding operation. With 
the above system the rate of data recording is lim- 
ited by the tape punch to about two 8-character data 
points per second. 

The possibility of erroneous recordings is greatly 
diminished by various self-checks incorporated into 
the system. The analog section of the system (cur- 
rent amplifier and voltage-to-frequency converter) is 
subject to the small errors inherent in analog devices 
and is therefore checked and calibrated several times 
during an experiment. The calibration is transferred 
directly to punched tape to preclude errors in record- 
ing. The nature of the digital section, and indeed of 
digital systems in general is such that the occurrence 
of small errors in transferring digits through the system 
is highly improbable. To ensure that such errors, if 
they occur, will be detected, an odd-bit parity check 
is made at the tape punch, and the recording process is 
automatically stopped if an added or dropped bit is 
detected on the tape. A second and independent 
parity check is made during the process of converting 
the tape to punched cards. 

2.4. Experimental Method 

A representative spectral line measurement con- 
sisted of the photoelectric recording of the arc inten- 
sity at each of perhaps 30 equally spaced positions 
across the width of the arc. This set of measurements 
was then repeated at each of usually 27 wavelengths 
spaced non uniformly across the spectral line — span- 
ning altogether about 10 times the maximum full 
halfwidth of the line. Finally, without turning off the 
electronics, recordings were made of the intensity of a 
calibrated tungsten strip pyrometer lamp at the same 
27 wavelengths. 

An intermittent arc traversing technique was used 
in which the motion of the arc carriage was momen- 
tarily stopped for an interval of time sufficient to 
overcome photomultiplier delay effects and to record 
the intensity at that point before resuming motion 
and proceeding to the next point. This eliminates 
the distortion of the arc profile which results from 
continuous movement of the arc carriage during the 
counting interval. Also, since the traverses always 
begin at precisely the same point (±0.001 cm), each 
recorded point corresponds to a given point in space 
referenced to a laboratory coordinate system, not a 
function of time or scanning rate. The recording 
positions are therefore repeatable and the recorded 
intensities directly comparable from traverse to 
traverse. 



The effect of integrating the intensity at a point for a 
finite interval of time is a temporal smoothing of noise 
at that point on the arc profile. Some further spacial 
smoothing occurs because of the necessity of observing 
a finite "slice" of the arc column. However, the ad- 
ditional loss of spacial resolution which would result 
from a continuous traversing or "scanning" of the 
arc image is effectively avoided by the intermittent 
traversing technique. In these experiments an inte- 
grating period of 0.1 sec was used. The width of the 
observed "slice" was determined by the slit width 
and in most experiments amounted to about 0.001 cm 
while the distance between samples was typically 
0.03 cm. 

The test gas — either pure argon or premixed mix- 
ture of argon and 1 percent hydrogen — was injected 
tangentially into the arc on both sides of the observa- 
tion chamber at a rate of 1.2 liter/min and a differential 
pressure of a few millimeters of water. Gas flow rates 
were monitored by flowmeters and adjusted to mini- 
mize flow in the observation chamber. Separate 
flows of pure argon at about 0.2 liter/min were injected 
tangentially at the electrodes. The difference in flow 
rates insured a net gas flow away from the observa- 
tion section, further precluding electrode vapor con- 
tamination by diffusion. 

The intensity measurements are calibrated with a 
tungsten strip lamp at a brightness temperature at 
6530 A of 2340.3 °K (IPTS 1948). The lamp was 
periodically calibrated by the Temperature Physics 
Section of the National Bureau of Standards. At 
the conclusion of each arc experiment it was placed 
at the arc position with the quartz window from the 
arc placed in the optic path to eliminate the necessity 
for determination of window transmission. A well- 
regulated d-c power supply maintained the strip lamp 
current at a constant level within 0.01 percent. The 
same punched program tape used to direct the mono- 
chromator through the sequence of wavelengths for 
the experiment was then used to record the calibrated 
strip lamp intensity at the same wavelengths as those 
of the experiment. 



3. Treatment of Data 
3.1. Abel Inversion 

In order to obtain the radial distribution of volume 
emission in the arc from side-on measurements it is 
necessary to solve the appropriate Abel integral 
equation [9] for each traverse across the arc. The 
technique employed is based upon one described by 
Freeman and Katz [10] in which the distribution of 
measured side-on intensities is fitted by a polynomial 
in the positional coordinate with origin at the center 
of symmetry. The solution of the integral equation is 
then obtained analytically for this polynomial and the 
volume emission at any desired point in the arc is 
calculated from this solution. By including odd powers 
in the fitting polynomial some minimal allowance is 
made for source asymmetry and for errors in locating 
the position of the arc "center of symmetry" [11]. 
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The polynomial fitting was accomplished by the 
use of the orthogonal properties of the Jacobi [12] 
polynomials Pty" a \xlR) whose inverse Abel transform 
is expressible in terms of P(?~ 1/2 ' 0) (2[r/#] 2 - 1) with odd 
n terms treated in the approximation of Freeman and 
Katz [11]. R is the distance from the arc center be- 
yond which the intensity is considered negligible. 
This technique results in the weighting of the data 
points by (1 — [x/R] 2 ) a . In the experiments reported 
here uniform weighting (a = 0) was arbitrarily used 
after experiments revealed no appreciable difference 
in the results obtained using a— (corresponding 
essentially to the Freeman and Katz method) and 
a= h (corresponding, except for the odd n terms, to a 
method described by Herlitz [13]). Individual arc 
traverses contained from 20 to 41 measurements at 
equally spaced x points and the resulting intensity 
profiles were fitted by the Jacobi polynomial expan- 
sions with 8 to 12 terms. 

3.2. Line Shape Parameters 

The volume emission intensities calculated from the 
Abel inversions were reassembled as a function of 
wavelength at each radial point and a suitable the- 
oretical spectral line shape was fitted by an iterative 
least squares technique. In the case of the argon lines 
a dispersion (Lorentz) shape superimposed upon a 
linear background was used as the fitting function: 



I = A+B(k-\o) + 



1 + [(X-X )/A] 2 



with A, B, C, \o, and A being adjusted to minimize 
the sum of the squares of the residuals. In the case 
of the hydrogen /§ lines the line shape fitting was pre- 
ceded by subtracting from each experimental Abel 
inverted wavelength scan the corresponding pure 
argon arc spectrum at the same radial point. This 
pure argon background spectrum was obtained as a 
part of each experiment in the same manner and under 
the same arc conditions as was the argon-hydrogen 
mixture data. The red and blue wings of the resulting 
difference spectrum were then averaged and this line 
shape was fitted to an appropriate theoretical H^ line 
shape: 

l = A+BS(a) 

where a= 1.25 X 10- 9 |A- \ \(N e )- 2l \ S(a) is a theo- 
retical line shape [14], and A, B, and the electron 
density A e (cm~ 3 ) are adjusted to minimize the residual 
function. The Griem, Kolb, and Shen tabulations of 
S(a) were fitted piecewise by cubics in all the tabulated 
intervals except the first where a parabola was used 

TO 

with the condition -^— = at a = Q. The asymptotic 
da 

wing formulas [15] were used beyond the range of 
the tables. In this manner piecewise analytic ap- 
proximations to S(a) were constructed for H/3. In 
addition to the dependence of the H# line shapes upon 
N e implied by the definition of a, there is a further 



weak dependence upon N e and T. This dependence 
was included in these calculations by the actual use 
of different piecewise line shapes calculated for 
N e = 10 17 and T= 13 000 °K, N e = 10 16 and T= 10 000 °K, 
and for N e = 10 15 and T=7500 °K-the first and last 
being obtained by interpolation in the Griem, Kolb, 
and Shen tables. 

All of the above calculations as well as division by 
the strip lamp calibration data and conversion to 
absolute intensities are programmed to be performed 
in one pass through the IBM 7094 digital computer. 
At each radial point the computer plots the experi- 
mental radially-resolved wavelength scans and tabu- 
lates all of the line fitting parameters; in particular 
the argon line centers A.o, 2 -half widths A, and total 
intensities 77 CA, the hydrogen line total intensity 
B JS(a)da, and the electron density N e . 

No slit width corrections were necessary for the 
hydrogen lines because of their relatively great width. 
However, for the argon lines slit or apparatus function 
corrections were required. For this purpose the pro- 
gram for dispersion line shape fitting was systemati- 
cally calibrated using various ratios of apparatus 
function width to dispersion line width and using 
various distributions of wavelength points of measure- 
ment across the line. The calibration was done for 
both Gaussian and triangular apparatus functions and 
resulted in the application of small slit width dependent 
corrections to the computed line widths and, for very 
wide slits, to the line intensities. This procedure 
should be quite reliable for lines which have nearly 
a dispersion shape. 

The apparatus functions observed with the 15 /x 
entrance and exit slits used in most of these experi- 
ments were approximately Gaussian with J -half widths 
(one-half of the full width at the half maximum inten- 
sity level) of 0.1 A. One experiment was performed 
with 150 fju slits in which the apparatus function was 
nearly triangular with a 2 -half width of 0.75 A. 

Examples of the line shape fitting are given in figures 
3 and 4. The AI 4300 A line can be seen to exhibit 
some skewness. This skewness together with the 
unsymmetric disposition of the wavelength points 2 
causes the dispersion shape fitting routine to produce 
a line shape which in this example is about 3 percent 
narrower than a good eye-fit would probably yield. 
The apparatus function for both lines in this example 
was the Gaussian shape with i -half width of 0.1 A. 
In figure 3 the computer-fitted o ^-halfwidth of 0.771 A 
required a correction of 0.013 A for apparatus function 
resulting in an estimate of 0.758 A for the ^-halfwidth 
of the spectral line itself. In figure 4 the intensity 
scale for H/3 is displaced vertically by an amount which 
has no particular physical significance being deter- 
mined by the subtracted hydrogen-free background 
spectrum. 

Beyond a radial distance of 0.2-.25 cm, the line 
shape fitting procedure was considered unreliable due 
to the sharpness of the lines and the attendent small 



2 The asymmetric distribution of wavelength points apparent in figure 3 is, to some extent, 
unavoidable unless an excessive number of uniformly spaced wavelengths is used because 
of the dependence of the line shift upon the radial position of the observation. 
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Figure 3. ^Experimental intensity points in the neighborhood of 
A I 4300 A after Abel inversion (0.13 cm from axis of 90 ampere 
are in argon with —7% H2) showing computer fitted background 
and dispersion shape: > -halfividth = 0.771 A. total line emission 
= 9.44 X 1 17 photons cm~ 3 sec~' 
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FIGURE 4. ^Experimental intensity points in the neighborhood of 
Hp 4861 A after Abel inversion {same conditions as fig. 3) showing 
computer fitted background level and theoretical line shape: 
electron density = 7.45 X 10 16 cm' 3 , total line emission = 8.10 X 10' s 
photons cm~ 3 sec' 1 . 



number of points defining their shapes. At these 
distances, the low number density of the emitting 
species causes a reduction in the signal/noise ratio 
with the effect that the determination of the true line 
shape from a statistically small sample of noisy data 
is improbable. For these reasons, the data from 
regions lying beyond an arbitrarily chosen radial dis- 
tance of 0.2 cm (0.16 cm for ionic lines) is not reported 
here. 

Argon line shifts and V2 -half widths corrected for the 
0.1 A Gaussian apparatus function are plotted in 
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Figure 6. Dependence of the red shift of the AI 4300 A line upon 
electron density. 



figures 5 to 10 against the electron density as deter- 
mined from the H^ line shapes. The wavelengths for 
the unshifted lines were obtained by extrapolation to 
N e = 0. Both the widths and shifts are seen to be 
approximately proportional to the electron density in 
rough agreement with theory [9|. The proportionality 
constants obtained by least squares fitting of this 
data are shown in table 1. 
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FIGURE 7. Dependence of \-halfwidth upon electron density for 
Al 6965/1 
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Figure 10. Dependence of the blue shift of the All 4806 A line 
upon electron density. 



Table 1. Stark 1/2 -half widths, w, 
and shifts, d, for argon lines 





w 


s [A cm 3 ] 


d 


[A cm 3 ] * 




Expt. 


Theor. 


Expt. 


Theor. 


A I 4300 A 
A I 6965 A 
All 4806 A 


0.100 
.048 
.0148 


0.124 
.061 
.0058 


0.103 

.039 

- .0073 


0.132 

.079 

(-.0022). 



*Shifts to the red are taken as positive. 



3.3. Transition Probabilities 

The calculation of the argon transition probabilities 
requires a knowledge of the plasma composition at 
each arc data point. For conditions of local thermo- 
dynamic equilibrium this can be computed in prin- 
ciple from the total pressure (1 atm), measured electron 
density, and total H^ line intensity assuming that only 
hydrogen and argon are present and using the theoret- 
ical H. transition probability, Dalton's law, the con- 
dition for quasi-neutrality and the appropariate mass 
action laws (Saha equations) for the possible ionization 
and dissociation reactions. The fact that the hydro- 
gen was present in trace amounts permitted some 
simplification of the calculations as has been dis- 
cussed previously [5, 16]. This simplified calculation 
requires a table of the composition of a pure argon 
plasma at atmospheric pressure and employs small 
corrections computed iteratively from the observed 
H^ line intensity. 

The pure argon plasma composition was machine 
computed by a method similar to that of Drellishak, 
Knopp, and Cambel [17] including first order Debye 
corrections. The partition functions as functions of 
the cutoff energy were analytically smoothed and made 
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to assume the asymptotic hydrogenlike form at the 
ionization energy limit (9). The number density of 
triply ionized argon atoms was assumed negligible 
thus limiting the calculations to temperatures below 
about 20 000 °K. The computed composition table 
including upper level populations pertinent to the 
present experiments is presented as table 2. Griem's 
treatment [18] of high density corrections to plasma 
composition calculations suggests that table 2 should 
be accurate to approximately 1 percent. 

Finally the electron density and Hp line intensity 
obtained from the H# line shape fitting are used itera- 
tively in conjunction with this table to determine the 
total hydrogen content and the number density of 
excited argon atoms at the radial points considered. 



The argon transition probabilities are then found by 
division of the total line intensities (in photons per 
cubic centimeter per second) by these number densi- 
ties at corresponding radii. The individual values 
of the transition probabilities obtained in this manner 
are presented in figures 11, 12, and 13. The averages 
of these experimental transition probabilities are 
given in table 3. 

The transition probabilities of the lines other than 
6965, 4300, and 4806 were determined relative to one 
or another of these three reference lines. The lines 
lying in the red and infrared were measured using 
the digital apparatus described above with full Abel 
inversion but with wider slits and with an abbreviated 
wavelength scan. The plasma composition was 



TABLE 2. Argon plasma composition at 1.00 atm pressure* 



T 


Z 


z, 


z 2 


Uo 


Ux 


V, 


Ne 


N 


V, 


N, 


5000 


1.0000 


5.3246 


7.904 


127106 


222812 


329954 


1.322E 


12 


1.468E 


18 


1.322E 


12 


3.642E-07 


5500 


1.0000 


5.3752 


8.034 


127101 


222803 


329940 


7.174E 


12 


1.335E 


18 


7.174E 


12 


1.441E-04 


6000 


1 .0000 


5.4188 


8.158 


127093 


222787 


329916 


2.944E 


13 


1.223E 


18 


2.944E 


13 


2.135E-02 


6500 


1.0000 


5.4567 


8.278 


127081 


222762 


329879 


9.746E 


13 


1.129E 


18 


9.746E 


13 


1.485E 00 


7000 


1.0000 


5.4901 


8.393 


127063 


222726 


329825 


2.724E 


14 


1.048E 


18 


2.724E 


14 


5.702E 01 


7500 


1.0000 


5.5196 


8.506 


127039 


222678 


329753 


6.648E 


14 


9.775E 


17 


6.648E 


14 


1.36 IE 03 


8000 


1.0000 


5.5459 


8.615 


127009 


222617 


329662 


1.453E 


15 


9.147E 


17 


1.453E 


15 


2.205E 04 


8500 


1.0000 


5.5695 


8.722 


126971 


222542 


329549 


2.898E 


15 


8.579E 


17 


2.898E 


15 


2.599E 05 


9000 


1.0001 


5.5908 


8.827 


126926 


222453 


329415 


5.352E 


15 


8.050E 


17 


5.352E 


15 


2.347E 06 


9500 


1.0001 


5.6101 


8.929 


126875 


222350 


329261 


9.254E 


15 


7.543E 


17 


9.254E 


15 


1.693E 07 


10000 


1.0003 


5.6276 


9.029 


126818 


222235 


329089 


1.51 IE 


16 


7.041E 


17 


1.51 IE 


16 


1.009E 08 


10500 


1.0005 


5.6437 


9.126 


126755 


222109 


328900 


2.342E 


16 


6.527E 


17 


2.342E 


16 


5.104E 08 


11000 


1.0009 


5.6584 


9.221 


126688 


221975 


328699 


3. 166 E 


16 


5.987E 


17 


3.466E 


16 


2.239E 09 


.11500 


1.0015 


5.6720 


9.314 


126618 


221836 


328491 


4.910E 


16 


5.41 2E 


17 


4.9 10E 


16 
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*Zo, Z\, and Z-> are internal partition functions for neutral, singly ionized, and doubly ionized argon atoms respectively. 
Uo, Ui, and U-> are effective ionization energies [cm -1 ]. 

N e , A 7 !), N\, and N 2 are number densities of electrons, neutral, singly ionized, and doubly ionized argon atoms [cm -3 ]. 
Columns headed by a wavelength contain number densities in the upper state of the line: e.g., under 4300 AI is tabu- 
lated N„g„ expi-Ejm/Zv with g„ = 5, £,, = 116999 cm" 1 . 
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FIGURE 11. AI 4300 A transition probability values. 

Ave. = 4.11 X 10 5 sec- 1 . Std. dev. = 0.29x 10 5 sec" 1 . 



FIGURE 12. AI 6965 A transition probability values. 

Ave. = 7.10xl0 6 sec" 1 . Std. dev. = 0.55x 10" sec" 1 . 



obtained from measurements of the AI 6965 A in- 
tensity using a transition probability of 7.10 X 10 6 
sec -1 . The lines in the blue were measured using 
an analog computer [19] for the direct recording of 
the Abel inverted spectra. In this way it was felt 
that the separation of overlapping lines (see fig. 1) 
could be better accomplished — although somewhat 
subjectively. The line intensities were determined 
by planimeter integration and the plasma composition 



was derived from the measured AI 4300 A intensity 
(using ,4 = 4.11 X 10 5 o sec _1 ) for the neutral lines and 
from the All 4806 A intensity (using ^ = 1.31X10 8 
sec -1 ) for the All 4348 A line. Because of the less 
satisfactory estimates of the line wing corrections, 
the limited number of experiments, and the uncer- 
tainty in the reference transition probabilities the A 
values for these lines are slightly less accurate than 
those for the 4300, 6965, and 4806 lines. 
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FIGURE 13. All 4806 A transition probability values. 

Ave. = 1.31 XlO 8 sec" 1 . Std. dev. = 0.21 X 10 8 sec." 1 . 
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Figure 14. The hydrogen/ argon atomic ratio, k, as a function of 
radius for some 60 ampere experiments. 
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3.4. Hydrogen- Argon Ratio 

The ratio of the total hydrogen content to the total 
argon content as determined from the H^ shape and 
intensity measurements alone is shown in figure 14 
as a function of radial position for some 60 A experi- 
ments. This ratio is defined as 

k = (2N H . 2 + N H + N H+ )l(N A + N A+ + N A+ + +. . .) 



where /V H2 , /V H , etc., are number densities of the species 
H 2 , H, etc. All of the experiments in figure 14 
utilized a premixed tank of the hydrogen-argon mix- 
ture with a value for this ratio of about 0.02 as deter- 
mined by a mass spectrometric analysis. This figure 
illustrates the "demixing" effect [9] which precludes 
the use of cold gas composition data in determining 
the plasma composition in arcs containing gas 
mixtures. 



4. Discussion of Errors 

The recent appearance of a detailed and self-con- 
sistent theoretical treatment [9] of the problems of 
laboratory plasma sources and of spectroscopic plasma 
diagnostics together with the relatively large data 
sample sizes that can be accumulated with automatic 
recording and computing techniques makes possible 
a much more thorough error analysis than has hereto- 
fore been attempted. A discussion of several sources 
of error follows. In general it is found that the ex- 
pected experimental and theoretical errors are of 
roughly comparable magnitudes. 

4.1. Electronic 

The maximum imprecision of the recorded intensi- 
ties due to instrumental uncertainties was determined 
to be of the order of 0.15 percent. This figure results 
from the treatment as independent errors of the maxi- 
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mum effects of instrument nonlinearities and zero 
drift over a period of 1 hr with the combined effects 
of ambient temperature variations of plus or minus 
5 °C and ± 10 percent line voltage fluctuations, and 
includes the contributions from the photomultiplier 
high-voltage power supply, the current amplifier, the 
voltage-to-frequency converter and the electronic 
counter. 

An additional variation has been observed in the 
photomultiplier dark current on the higher gain range 
settings of the amplifier. This effect became notice- 
able after the dark current had been suppressed and 
appeared to consist of a long-term drift of the dark 
current of as much as 1 percent of full scale accom- 
panied by short-term random noise of a similar mag- 
nitude. 

The effects of photomultiplier fatigue on the 
photocurrent measurements were made negligible 
by restricting anode currents in the photomultiplier 
to less than 10" 7 A, a value which previous investiga- 
tions had shown to be below the threshold of fatigue 
effects for these tubes. 

A certain amount of delay in the response of the 
photomultipliers to step inputs was observed. A 
rise time of about 0.2 sec was measured for the photo- 
tubes with S-20 response, while the delay was un- 
measureable for the S-4 phototubes, implying some 
dependence on photocathode material. To ensure 
sufficient time for any effects of photomultiplier delay, 
mechanical relaxation or transient vibrations to decay, 
the recording of the photocurrent was postponed until 
0.25 seconds after the arc chamber had completed 
its shift to a new position. 



4.2. Calibration 

The absolute calibration of the source intensities 
with a tungsten strip lamp is subject to a maximum 
uncertainty of about ± 5 percent. The brightness 
temperature calibration of the strip lamp at NBS is 
qualified by a certificated uncertainty of ± 7 °C 
relative to the IPTS which would produce uncertain- 
ties of about ± 3 percent in the transition probabilities 
of both neutral and ionic lines. For the calculation 
of the strip lamp emission at wavelengths other than 
the calibration wavelength (6530 A) the tungsten 
emissivities given by DeVos [20] were used. These 
are estimated to be accurate to within ± 2 percent. 
The maximum uncertainty due to this cause would 
therefore presumably be ± 2 percent in the blue and in- 
frared spectrum and would decrease to zero at the 
calibration wavelength. 

The errors in intensity due to uncertainties in 
determination of the correct strip lamp current are 
around ±0.07 percent, almost two orders of magnitude 
below the errors discussed above and therefore 
negligible. The combined expected uncertainties are 
therefore ±5 percent for blue and infrared lines 
decreasing to ±3 percent for lines in the red. 



4.3. Arc Source 

The argon used in these experiments was technical 
grade argon with typically a purity of 99.996 percent 
containing mainly nitrogen and oxygen as contami- 
nants. The gas flow system in the arc chamber is 
designed to prevent atmospheric and electrode con- 
taminants from entering the observation chamber and 
indeed no impurities were observed spectroscopically 
during the experiments. From tins observation the 
impurity upper limit and corresponding transition 
probability error of the order of 1 percent can be set 
based upon a minimum detectable signal and estimated 
transition probabilities for lines of the likely impuri- 
ties, C, N, and 0. 

The arc current regulator maintains the arc current 
constant to within 0.05 percent for periods of 1/4 or 
1/2 hr. Effects of a long-term drift presumably due 
to ohmic heating of the oil-cooled regulator shunt 
were minimized by monitoring the arc current at 
intervals with a potentiometer and precision shunt 
and making small manual corrections. Oscilloscopic 
observation of the arc current reveals a 0.3 percent 
360 cycle ripple when the rotary converter power 
supply is used. Assuming the totally unrealistic 
worst possible case — that the averaging over the 0.1 
sec counting time results in a measurement of the elec- 
tron density at one current extremum and of the line 
intensity at the other — then the ensuing error in the 
transition probabilities would be of the order of 0.5 
percent and 1 percent for the AI and All lines respec- 
tively. 

The arc positioning mechanism is capable of locating 
the arc chamber to within 0.002 cm, with a repeat- 
ability an order of magnitude better. Since these 
positioning errors are attributed to the lead screw 
and are therefore noncumulative, they will appear 
on the recorded profile of the traverse as noise amount- 
ing at most to about 0.4 percent rms of the maximum 
signal. 

The assumed total arc pressure of 1 atm is clearly 
an approximation. No record was kept of the prevail- 
ing atmospheric pressure during the experiments so 
that day-to-day pressure fluctuations of 2 or 3 percent 
must be expected. The over-pressure of the argon 
in the arc required to maintain the steady gas flow 
was 0.0003 atm. The arc column over-pressures due 
to magnetic self-compression and due to the mo- 
mentum transport accompanying the heat transport 
can be estimated [9] for the 90 A arc as 0.0003 atm 
and 0.02 atm respectively. Thus the actual pressure 
at all points in the arc was at all times probably within 
5 percent of the 1 atm assumed for the calculations. 
By repeating the plasma composition calculations 
(table 2) at other pressures the effect of this pressure 
uncertainty upon the transition probabilities can be 
evaluated. The result is that a ±5 percent error in 
pressure produces a+ 1 percent error in the neutral 
line transition probabilities and a ±8 percent error 
in the ionic line values. 
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4.4. Optical 

Any nonuniformity of the arc in the axial direction 
within the observation chamber will result in transi- 
tion probability errors due to the finite height of the 
spectrometer slit which collects light from a region 
0.05 cm long in the axial direction and because of the 
resulting signal averaging which will be different for 
line widths, neutral line intensities, and ion line in- 
tensities. Estimates of the maximum errors from this 
source based upon measurements taken during a 
crude scan along the arc axis are 0.5 percent and 
4 percent for neutral and ionic line transition prob- 
abilities respectively. Similar errors arise in the steep 
gradient regions of the arc because of imaging errors 
of the optical system and the finite width of the 
entrance slit. The combined lateral image smearing 
was determined to be less than 0.005 cm for the 15 /x 
slit by moving the edge of a tungsten strip filament 
lamp across the optic axis at the arc position and ob- 
serving the "halfwidth" of the signal step. The 
corresponding maximum transition probability errors 
in this case are also 0,5 percent and 4 percent respec- 
tively at the edges of the arc. In all cases the experi- 
mental values should be too high. Turbulence, if 
present in the arc, would lead to similar errors. How- 
ever, estimates of the Reynolds number for the condi- 
tions of these experiments lie well within the laminar 
flow region suggesting that this should not be a 
problem. Optical refractivity gradients in the arc 
are expected to cause a negligible distortion of the arc 
intensity profile even in the vicinity of strong lines 
[9]. This expectation was partially confirmed experi- 
mentally by the failure to detect a significant displace- 
ment of the apparent position of the edge of a tungsten 
strip lamp filament imaged into the arc plasma at any 
observing wavelength near the AI 6965 A line. 

Scattering of light by the optical elements and arc 
chamber windows would cause an increase in the 
apparent transition probabilities at the higher radii, 
i.e., at the edges of the luminous arc column. Since 
a slight variation of this type was indeed observed, 
several experiments were performed with no window 
in the arc chamber. These failed to show a signi- 
ficant difference at the arc edges from those in which 
the window was included in the optical path. Direct 
rough measurements of the light scattered by the 
optical elements showed this amount to be less than 
0.3 percent of the maximum source intensity, which 
could not possibly account for the observed radial 
variation of the ^-values. 

The shot noise arising from statistical fluctuation 
in the arrival of photons at the photomultiplier can 
be estimated from the source intensity and the geo- 
metry of the optical train. It is estimated to be always 
less than 0.5 percent for the strip lamp calibration 
and less than 0.1 percent for the arc, except at the 
extreme edges. 

For some experiments the limiting aperture of the 
optica] system was reduced to one-half its original 
size. Any observed difference in the results of these 
experiments would have pointed out errors and 



misalinements of the optical system. No such 
differences were observed, indicating that the optical 
system was properly alined and that the original 
aperture was indeed limiting. 

The degree of self-absorption can be estimated by 



calculating the optical depth t = 2 \ k' 

Jo 



and 



K' = e/B(T) 

B(T) = 2hc 2 k-He hc ' XkT -iy i 



k' is the effective absorption coefficient (cm -1 ), and e 
is the observed radially resolved emission coefficient 
per steradian (erg sec -1 cm" 3 sr _1 cm -1 ). The results 
of such calculations at the centers of the three lines 
at 4300, 6965, and 4806 A are respectively 0.002, 0.11, 
and 0.002 for the optical depths, r, in the worst cases 
encountered. Since the Abel inversion is most un- 
stable near the center due to the infinite kernel at that 
point it is to be expected that the consequences of 
errors in the side-on intensity measurements due to 
self-absorption will normally be most severe in the 
radially inverted values near the center of the arc. 
This is borne out by exact analytic calculations of 
simple models and by the approximate Freeman and 
Katz [10] treatment of more realistic models. These 
calculations indicate that for small degrees of self- 
absorption and for bell-shaped profiles the relative 
error in the radially resolved emission coefficient is of 
the order of KqR at the arc center where k' () is the effec- 
tive absorption coefficient at the center and R is, as 
before, the nominal arc radius. Thus only in the case 
of the AI 6965 A line will self-absorption play a role 
leading to transition probabilities that are too low 
by about 5 percent (k' () R averaged over the line shape) at 
the arc center and increasing to the correct value at 
the edge. To some extent this conclusion appears 
to be supported by o the data of figures 11 and 12 which 
show for the 6965 A line a more pronounced depression 
of the transition probabilities at high electron densi- 
ties near the arc center. 

4.5. Analytical and Theoretical 

In general, the accuracy of an Abel inversion based 
on curve fitting techniques is dependent upon the form 
of the intensity profile of the arc traverse and the fitting 
function, the number of polynomials used in the fitting 
function, the number of increments examined during 
each arc traverse, and the noise present on the in- 
tensity profile. To assess the effects of these variables 
on the accuracy of the inversion, analytical test func- 
tions were constructed with varying amounts of super- 
imposed noise drawn from a table of random normal 
deviates. These test functions, which are of the same 
general shape as the observed profile of the arc tra- 
verse, were then inverted by the machine program used 
in processing the experimental results. As expected, 
for data with a relatively high noise level a small num- 
ber of polynomials in the fitting function produces the 
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best inversion in terms of accuracy of fit and minimum 
noise amplification. Conversely, no significant im- 
provement in accuracy results from using large num- 
bers of increments and terms unless the signal/noise 
ratio of the data is quite high. The accuracy of the 
inversion is, however, relatively insensitive to varia- 
tions in number of increments for noise levels greater 
than 0.3 percent due to the smoothing inherent in the 
curve fitting process and in this range the relative 
noise level of the inverted profile is almost identical 
to that on the original recorded intensity profile. 

The results of tests on analytical functions were 
borne out by further tests of the inversion on actual 
experimental profiles. The number of data points of 
representative many-point arc traverses was halved 
and quartered by systematically discarding alternate 
data points, and the resulting "thinned" traverses 
inverted and compared to the originals. These tests 
also showed that, at least for the data encountered 
here, no further accuracy could be attained by taking 
more than 20 or 30 measurements per arc traverse. 
Similarly, the use of more than 8 or 10 terms in the 
polynomial expansion was again not warranted and, 
in fact, increased the noise level in the results. 

Numerical tests of the dispersion shape fitting pro- 
gram carried out by superimposing noise upon calcu- 
lated dispersion shapes reveals that the relative rms 
noise on the computed line intensity and halfwidth 
is roughly equal to the relative noise in the input 
intensity data at the half maximum intensity level. 
Thus, little smoothing of the noise is accomplished in 
this step. This is due mainly to the influence of the 
noise upon the background or continuum level de- 
termination by the program. In these experiments 
this determination was dependent upon a compara- 
tively small effective number of points due to the 
wavelength point distribution employed. 

The use of the area of the fitted dispersion shape 
as the spectral line intensity obviates the necessity 
for correcting measured line intensities for that portion 
of the line lying below and beyond the end points of 
the wavelength scan. This is a significant correction 
amounting to 13 percent for scans of only ten half- 
widths and may account for part of the difference 
between the transition probabilities observed here 
and the lower values of other workers. Since the 
experimental line shapes even after correction for 
apparatus functions are not expected to be quite 
dispersion shapes [9] it was necessary to ascertain 
the effect of departures from the dispersion shape 
upon the linewidth, shift, and intensity values as meas- 
ured in these experiments. For this purpose the 
dispersion shape fitting program was applied to a 
theoretical complete asymmetric line shape [21] which 
appears to be reasonably representative of the ex- 
perimental line shapes encountered here. The com- 
puted intensities and halfwidths were consistently 
low by 3 percent and 2 percent respectively. The 
line shifts were consistently overestimated by a mag- 
nitude approximately equal to the error in the *- 
halfwidth. No attempt was made to calibrate the 
dispersion shape fitting program for these theoretical 



line shapes because of the unavailability of suffi- 
ciently extensive tables of the necessary line shapes. 
Consequently, systematic errors can be expected in 
the results amounting to perhaps 5 percent and prob- 
ably always in the direction indicated by this example. 

At the heart of these measurements is the deter- 
mination of the electron density by the FL line shape 
measurement. Unfortunately a significant uncer- 
tainty in the results is associated with this step. 
For example, Griem [9] estimates that the electron 
densities obtained from FL line shapes are probably 
reliable only to about 5 percent based upon an anal- 
ysis of the theoretical approximations. In addition, 
as indicated in figure 4, some error is involved in the 
fitting of a theoretical line shape to the experimental 
data. Typically the rms deviation of the experimental 
points from the fitted curve is 2 or 3 percent of the 
maximum H^ intensity. (This is somewhat better 
than the dispersion shape fitting to the argon lines 
which typically has an rms deviation of 3 to 5 percent. 
However, in part this greater error is explained by 
the influence of the apparatus or slit function which 
is not taken into account in the fitting and by the ap- 
proximation of the true theoretical shapes by the dis- 
persion shape.) The subtraction of the pure argon 
spectrum from the argon + hydrogen spectrum re- 
moves the major part of any uncertainties in the FL 
shape measurement due to irregularities in the argon 
continuum underlying the hydrogen line. Any errors 
in this step due to small changes in the radial tempera- 
ture profile caused by the addition of the trace of 
hydrogen are expected to be negligible but pre- 
sumably are included in this 2 or 3 percent deviation. 
As a partial check of these FL measurements in two 
of the 60 A experiments H a line shapes were also 
recorded. From the electron density, temperature, 
and hydrogen to argon ratio as determined from the 
H^ measurements the theoretical intensity and shape 
ol the H a lines were computed using the Griem, Kolb, 
and Shen tabulations. Figure 15 shows a compar- 
ison of an experimental H a line shape (after the Abel 
inversion and subtraction of the pure argon spectrum) 
with the theoretical shape based upon the FL meas- 
urements at the same radial point. The agreement 
is seen to be satisfactory in the wings and at the center 
but the halfwidth of the theoretical shape is 10 percent 
or so smaller than the experimental one. This dif- 
ference appears to be in accord both in magnitude 
and direction with the correction to the theoretical 
H a shapes which Griem et al., [14] estimate would 
accrue by taking into account the improved ion field 
strength distribution functions of Mozer and Baranger. 
Thus, in conclusion, it is reasonable to assume that 
the electron densities obtained in these experiments 
from the H^ line shapes should not be in error by much 
more than about 5 percent. This results in 8 and 20 
percent uncertainties respectively in the AI and All 
transition probabilities. 

Griem [22] has discussed in detail the validity of the 
local thermodynamic equilibrium (L.T.E.) assumption 
in laboratory plasmas. Using his criteria one finds 
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FIGURE 15 . Experimental intensity points in the neighborhood of 
H tt 6563 A after Abel inversion (on axis of 60 ampere arc in argon 
with ~ 1% H 2 ) with theoretical H« line shape computed from 
parameters obtained by H# measurements: N e = 7.35X 10 16 cm~ : \ 
k = 0.0197. 



that for homogeneous plasmas of the dimensions of 
the arcs of these experiments complete L.T.E. would 
be probable down to the lowest electron densities 
actually utilized in the measurements — 2 X 10 16 for 
neutral species and 5 X 10 16 for ionic species. Again 
following Griem's treatment one finds that the greatest 
average distance traveled by any species in the arc 
between equilibrating collisions is that of the ground 
state neutrals of the order of 10~ 2 cm at the axis. 
Thus, the measured electron density and line intensi- 
ties at a particular point in the arc are not to be inter- 
preted as corresponding exactly to L.T.E. conditions 
at the local electron temperature at that point but 
rather as appropriate to some kind of averaging of 
near-equilibrium distributions which exist over a 
distance of the order of 10 -2 cm. The use of this 
distance, then, permits some estimate of the effects 
of the temperature and composition gradients present 
in the arc upon the accuracy of the transition proba- 
bility determinations. Clearly an over-estimate of 
the errors is obtained by coupling line intensity meas- 
urements at one point with equilibrium calculations 
based upon electron density measurements 10~ 2 cm 
away. Such a procedure results in error estimates of 
1 percent for the neutral lines and 2 percent for the 
ionic lines at the center of the 90 A arc. Away from 
the center this equilibrating collision distance is not a 
proper measure of the L.T.E. smearing because of the 
additional equilibrating influences of the diffusion of 
ions and excited atoms from the hot core and the ab- 
sorption of resonance radiation from the core. The 
fact that the measured transition probabilities show no 
strong dependence upon the electron density, however, 
suggests that if departures from L.T.E. are not serious 
at the arc center they are probably also not serious at 
least out to 0.2 cm where the lowest electron density 
measurements were taken. 



5. Discussion of Results 

The estimated transition probability uncertainties 
resulting from this error analysis assuming the errors 
are independent are included in table 3 together with 
transition probability values obtained by other investi- 
gators. Most of the errors which have been discussed 
are systematic. However, a few — the electronic 
noise, the numerical curve fitting errors, the day-to- 
day pressure variation, and the dark current drift 
errors — are not and will tend to be reduced in magni- 
tude by the averaging over the many experiments or 
radii. In arriving at the uncertainty estimates in 
table 3 due allowance for this kind of error was taken, 
where appropriate, by division by the square root of 
the sample number. A similar calculation applied to 
the neutral line width and shift measurements gives 
an estimated 6 percent uncertainty. 

The three dominating errors in the average transition 
probability values determined here are the absolute 
intensity calibration errors, the systematic errors in 
the argon line shape fitting due to apparatus functions 
and asymmetry, and the estimated uncertainty in the 
hydrogen line broadening theory. It is now possible 
to reduce the effect of calibration uncertainties from 
the 3 to 5 percent estimated in these experiments 
to about 1 to 2 percent, by means of new NBS radio- 
metric calibration services which have become avail- 
able on a limited basis [24] since these experiments 
were undertaken. Also, with some effort, the argon 
line shape fitting can be improved by the calculation 
and use of more extensive theoretical line shape tables 
including asymmetry and possibly with apparatus and 
Doppler broadening properly folded in. Thus the 
experimental errors can probably be reduced by a 
factor of two without excessive effort. However, 
without a concomitant improvement in the accuracy 
of the electron density determination these experi- 
mental improvements would only marginally reduce 
the estimated overall errors. The uncertainty in this 
step is the major contribution to the error estimates 
in table 3. Griem's estimate of the accuracy of the 
hydrogen line shape calculation is not significantly 
changed by consideration of the small amount of rela- 
tively contradictory experimental evidence [9|. Fur- 
ther experimental or theoretical studies should permit 
a better assessment of the accuracy of this extremely 
useful plasma diagnostic technique. 

The origin of the apparent decrease of the transition 
probabilities with increasing electron density (see 
e.g., fig. 11) is unknown. Of the error sources de- 
scribed above, the only pertinent ones estimated to 
be of sufficient magnitude to encompass this effect are 
hydrogen line broadening theory errors, line shape 
fitting errors, and possibly departures from L.T.E. 
This deviation from a constant transition probability 
agrees qualitatively with the discrepancies observed 
by Wiese et al. [25] in pure hydrogen arcs. 

The scatter of the transition probability values (figs. 
4 to 13) is in agreement with the noise level estimates 
described above. These were mainly electronic and 
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amounted to 1 or 2 percent of the maximum signals. 
When expressed in terms of the half -intensity of a line 
this relative noise level is higher by a factor of two or 
more, and, consequently, since neither the Abel inver- 
sion nor the dispersion line fitting seems to modify the 
noise level, the line intensities should exhibit a scatter 
of several percent. When combined with a similar 
noise in the electron density values a resultant scatter 
of the transition probability values of the order of 5 to 
10 percent is to be expected. 

The theoretical transition probabilities of Garstang 
[23], the "hook" method values of Pery-Thorne and 
Chamberlain [1], and the arc measurements of Drawin 
[4] appear to agree generally with the present values in 
the sense that the error ranges assigned by the investi- 
gators overlap the present ones. (The theoretical 
values are good to perhaps a factor of two, Pery- 
Thorne and Chamberlain estimate an accuracy of 20 
to 30 percent for their values, and Drawin assigns 
error limits ranging from o 10 percent for AI 4044 A to 
20 percent for AI 4300 A.) Drawin's measurements 
were carried out in an arc containing a hydrogen-argon 
mixture. However, since his diagnostics were based 
upon hydrogen line intensity measurements and the 
use of the cold gas composition there is a possibility 
of a systematic error due to demixing. 

Gericke's argon transition probabilities [3] are based 
upon an argon arc temperature profile derived from 
measurements of the All 4348 A line intensity using a 
transition probability value of 68 X 10 6 sec -1 . The 
value obtained in the present work for this transition 
probability is 240 X 10 6 sec -1 . If this is correct then 
Gericke's transition probability values should be in- 
creased by about 35 percent 3 and would agree well with 
the present values. These same comments apply to 
Richter's value for the AI 4044 A line (27) which was 
evidently obtained similarly. The corresponding 
adjustments to Gericke's electron density values result, 
in the case of the AI 4300 A line, in values (in A) for 
the ^-halfwidth and shift both approximated by 0.12 
X 10 -16 [W e ] where N e is the number of electrons per 
cm 3 . This compares reasonably well with the present 
values (table 1) of 0.10 X 10 -16 [A^]. Thus the measure- 
ments reported here are evidently well in accord with 
Gericke's on a relative basis but the absolute results 
differ through the use of the H# line shape for the arc 
diagnostics rather than the All 4348 A transition 
probability used by Gericke. 

In the case of Olsen's results an estimate of the 
degree of self-absorption at the center of the All 
4806 A line 4 under his arc conditions leads to values 
for the optical depth, r, of 0.7 or 0.3 depending upon 
whether the line halfwidth is estimated from Griem's 
theoretical calculations [9] or is extrapolated from the 
experimental data of figure 9. In either case the 
magnitude is sufficient to cause a serious error in the 
radially resolved emission coefficient and to raise 
doubts about whether Olsen's observed apparent 



off-axis maximum was not in fact a consequence of 
this absorption rather than a manifestation of the 
Larenz and Bartels [28] "normal" temperature effect. 
In this event Olsen's temperature determinations and 
transition probabilities are open to radical reinter- 
pretation. 

Table 1 compares the measured values of the line 
widths and shifts (averaged assuming they are pro- 
portional to the electron density) with theoretical 
values computed from Griem's tables [9] for a tem- 
perature of 12 000 °K. Fortunately the theoretical 
values of the reduced widths and shifts are not strongly 
temperature dependent so that the difference between 
12 000° and whatever effective temperature is intro- 
duced by the averaging of the experimental values will 
cause a difference of only a few percent in the tabular 
entries. For the neutral argon lines the agreement 
of theory and experiment is roughly within the ex- 
pected accuracy of 8 percent for the experiments and 
20 percent for the theory. The ion line results, how- 
ever, are in very poor agreement with theory. For 
the line shift this is understandable since the theo- 
retical Stark shift for the line does not include the 
plasma polarization blue shift [9] which rough calcu- 
lations indicate will be larger than the Stark shift. 
For the width of the All 4806 A line errors due to both 
Doppler broadening and the slit width correction pro- 
cedure are important but do not appear to be sufficient 
to resolve the discrepancy between theory and experi- 
ment. The consequences of neglecting the Doppler 
broadening were estimated by folding the Doppler 
line shape (Gaussian with ^ -halfwidth of approxi- 
mately 0.03 A) into the apparatus function to give a 
broadened effective apparatus function. The use of 
this modified apparatus function (which significantly 
affects only the narrow ionic line) reduces the widths 
of figure 9 by approximately 0.005 A and decreases 
the proportionality constant in table 1 by about 5 
percent. The error in this proportionality constant 
stemming from the application of the slit width cor- 
rection procedure to the All 4806 A line where the 
line and apparatus function widths are comparable 
(line width: 0.08 A to 0.15 A: apparatus width: 0.1 A) 
is estimated at 10 or 20 percent. This figure is based 
upon the results of the 150 /x slit width experiment 
in which jthe apparatus function was wider than the 
AI 4300 A line at all radii. In this case application 
of the slit width correction procedure to the com- 
puter fitted line widths resulted in an average value for 



N, 



X 10 16 of 0.089 compared to the average 0.100 for 



3 This includes a correction discussed by Griem (26) to the ionization potential lowering. 

4 There appears to be an error in the equation for self-absorption in reference 2 (eq (6)). 
' n should probably be replaced by \ tt. 



all the narrow slit experiments. 
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